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Kondo effect in Dirac systems 
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Tsukuba Central 2, 1-1-1 Umezono, Tsukuba 305-8568, Japan 

We investigate the Kondo effect in Dirac systems, where Dirac electrons interact with the localized 
spin via the s-d exchange coupling. The Dirac electron in solid state has the linear dispersion and 
is described typically by the Hamiltonian such as Hk = nk ■ a for the wave number k where cr, 
are Pauli matrices. We derived the formula of the Kondo temperature Tk by means of the Green’s 
function theory for small J. The Tk is determined from a singularity of Green’s functions in the form 
Tk — D exp(—const. / p\J\) when the exchange coupling |J| is small where D = D/yjl + D 2 /(2p) 2 
for a cutoff D and p is the density of states at the Fermi surface. When \p\ <C D, Tk is proportional 
to |/r|: Tk — \p\ exp(—const./p|J|). The Kondo screening will, however, disappear when the Fermi 
surface shrinks to a point called the Dirac point, that is, Tk vanishes when the chemical potential 
p is just at the Dirac point. The resistivity and the specific heat exhibit a log-T singularity in the 
range Tk < T \p\/kB- Instead, for T ~ 0(|/r|) or T > |/r|, they never show log-T. 


I. INTRODUCTION 

Recently, the Dirac electron in solid state has been in¬ 
vestigated intensively because the Dirac cones in metals 
or semimetals have been realized in solids [l|-|hj. In sev¬ 
eral materials such as Graphene [bI4Io| as well as Bismuth 
compounds, there appear the conduction bands with the 
linear dispersion that are described by Dirac Hamilto¬ 
nian. The Dirac cone in the band structure appears in 
the surface states of topological insulators [ill 03 ■ Three- 
dimensional Dirac semimetals have also been realized in 
Na 3 Bi[Hj] and CdqAs9fl4l|. 

The Kondo effect, that occurs as the result of the ex¬ 
change interaction between dilute magnetic impurities 
and conduction electrons, is one of the most important 
phenomena in solid state physics [ill [HI- It is interest¬ 
ing to examine how the Kondo effect emerges in Dirac 
systems. The linear dispersion of Dirac electrons would 
affect the characteristic features of the Kondo effect. It is 
not even trivial whether the Kondo effect indeed appears 
in Dirac systems. The purpose of this work is to show 
unique and interesting properties of the Kondo effect in 
these systems. 

The s-d model in Dirac systems is closely related to the 
pseudogap model of the Kondo problem where magnetic 
impurities couple to conductive fermions with a density 
of states p(w) oc \uj\ r (0 < r) [Ml. In the pseudogap 
model, the density of states vanishes at the Fermi level. 
This feature is common to the Kondo problem in Dirac 
systems when the Fermi surface is point like, namely, the 
chemical potential is just at the Fermi point p = 0. It has 
been suggested that there is a phase transition between 
the local-moment phase and the strong-coupling singlet 
phase in the pseudogap Kondo and Anderson models [18[. 
A Dirac system is also interesting from the viewpoint of 
topology because the index theorem has been proved for 

Dirac operators (13, (HI 

In this paper we investigate the s-d Hamiltonian in a 
Dirac system by means of the Green’s function theory 
and evaluate the Kondo temperature Tk- We employ 


the decoupling procedure to obtain a closed solution for 
a set of equations for Green’s functions. Although the 
decoupling procedure is valid only for small J in the re¬ 
gion, T > Tk, this method is useful to derive the Kondo 
temperature. It turns out that Tk is crucially dependent 
on the chemical potential p. Tk vanishes for p = 0 at 
the Dirac point in our method. This shows the absence 
of Kondo screenings at p = 0 and suggests that there 
may be a transition at some critical value of J = J c 
from the local-moment phase to the singlet-formation 
phase at p = 0. In Dirac systems Tk is finite and is 
given by the formula being similar to that in the con¬ 
ventional Kondo effect: Tk oc D exp(—const./pi J|) when 
p|J| is small where p is the density of states and J is 
the exchange coupling constant. When p|J| is large of 
order 1, Tk is given by an algebraic function of p\J\ like 
Tk oc |/z|Q(p|Jj) for a some function Q(x). 

The paper is organized as follows. In section II we 
show the Hamiltonian. In section III Green’s functions 
are defined and a set of equations for them are derived. 
In section IV the Kondo temperature is derived. In sec¬ 
tions V and VI, we evaluate the electric resistivity and 
specific heat, respectively, and discuss their log-T singu¬ 
larity. The last sect ion is devoted to a summary. 

II. HAMILTONIAN 

We consider a massless Dirac Hamiltonian (Weyl 
Hamiltonian): 

Hd = ^2^1 ( V X k X (J X + VykyGy + V Z k z G Z ) , (l) 

k 

where 

*-(£) (2) 

with the annihilation and creation operators Cko- and , 
respectively. Here v x , v y and v z are velocities of conduc¬ 
tion electrons. a x , a y and a z are Pauli matrices. In our 
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model the Dirac electrons, described by this Hamiltonian, 
interact with the localized spin. The total Hamiltonian 
is written as H = H$ + H s d where 

H 0 = 'Y' J [{vxk x - iv y ky)c\,^c 14 + ( v x k x + iv y k y )c^c k f 

k 

+ ^/Cz( 4 t c kt - c^Ckj.) - /x( 4 T c kf + 4 c ki )], 

(3) 

Had. = ~ 2 lv E^ z ( C kt Ck 't _ C k 4 . Ck 'd) + ‘ 5 , + c ik 4 . c k't 

kk' 

+ S , _4 t c k / 4 .]. (4) 


N is the number of sites and we have included the chem¬ 
ical potential /x. S +, S_ and S z denote the operators of 
the localized spin. The term H s d indicates the s-d inter¬ 
action between the conduction electrons and the localized 
spin, with the coupling constant ■/EMI- J is negative, 
as adopted in this paper, for the antiferromagnetic inter¬ 
action. The Dirac Hamiltonian resembles the s-d model 
with the spin-orbit coupling of Rashba type 1221. We use 
the Green’s function method, following Ref.|22| to eval¬ 
uate the Hondo temperature and its related properties. 


III. GREEN’S FUNCTIONS 

We define thermal Green’s functions j23l. [2~il : 



G kk ' CT (r) = 

-{T T Cka(r)cl ltJ (0)), 

(5) 


H kk /(r) = 

— (T T Ckl ( T ) c k't (0) ), 

(6) 


5 ckt; c k't)) r = 

— (TtS z Ckt ( r )4't (0)) , 

(7) 

((S- 

-Ck|! C k /f))r = 

~(T t S- c k j_(r)4/^(0)), 

(8) 

«£ 

2 ckt; c k't)) r = 

C k -|-(T) C jj.,^ (0)) , 

(9) 

((S- 

-Ckt! c k /^))r = 

— {T t S- c k j,(r)4/^(0)). 

(10) 


Here T t is the time ordering operator. The Fourier trans¬ 
forms are defined as usual: 

GWo-fj) = — ^ e l “ nT G kk / CT (tw n ), (11) 

" n 

f™{t) = ^j2 e ~ iulnTF ™( iuj ^’ ( 12 ) 

" n 

((ffzCkt; c k't))r = ~n El e lUlnT c\ L ,^))i uln , (13) 

“ n 

((S.c^cl^r = W e ~ iUnT «S-C^<i',))i^ (14) 

" n 

We also use the following Green’s function, 

r k k'f(r) = -g E e~ laJn T kk /^(za; ra ) 

" n 

= ((5'zCkt;4, t )) T + ((-S_c k ^; 4't)) r ' (1^) 


We start from the commutation relations: 


[Ho, c kt ] 

[^o,c k 4_] 
[H s d, Ckf] 

[H s d, c k ;] 


-{v x k x ~ ivyky)cki + ( -v z k z + /x )c kt , 


(16) 

~{v x k x + iv y ky)c kt + (v~k z + /x)c k 4., (17) 

-^E(- 5zCk 't-^ Ck 't) 5 (18) 

k' 

(1 9 ) 

k' 


The equations of motion for G kk "|-(r) and F kk / are 


^G kk / t (r) = -5(r)<5 kk / + (~v z k z + /x)G k k't(r) 

- iy x k x - iv y k y )F kk >(T) 

+ 9^E[(( 5zCi 5t; c k't))^ + « s, - c 94.;4't»4 

q 

( 20 ) 


^H kk /(r) = (v z k z + /i)F kk / - (v x k x + iVykyjGw^r) 

~ 2N EE^Gif; c k't^ T ~ ((<S , +Cqt; 4't))r]- 

9 

( 21 ) 


Then the equation for G kk ^ reads 

{lOJn + /X - V z k z )Gw^(iu n ) = £) kk ' 

+ iy x k x tfy ky )7 kk * (ioJ n ) ^ y ^ ) rq k , 'f'(^n). (22) 

q 


The equations for ((S^Ckt; 4't)) an< l (('S’-CkU 4q-)) are 
similarly obtained as 


V z k z T /x) ( (S z C k ^, C^,^) )iuj n {^z^kk' 

T (v x k x iv y ky){{S z c^-, c k /^)) iu) n 

~ 2 TV El [((‘-’z c qt; 4't))®“n 4” 2 ((‘8 , ~ Cq 4’ c k"|-)) i “n 

q 

_ 2iV E [((‘8 , + Ck t C qJ. C q't: 4'f)) iaj ™ 

qq' 



(23) 
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(iid-n H - Vzkz H“ AO ((*^— OcE ^k / 'f‘))i w n 
— ( Vxk'x H - iVyky )(($_ Ck'i'^ ^-'k / 'i'))* Ci; n 

~ 2 ~nY 1 [((‘ S +‘ S '- c qt; c k't))“« + ^(( s - c q4-; c k't))^ 

q 

~ 2 iV [^‘^- Ck -l- C qt C q't’ C k't)) iw n 
qq' 

~ ((‘^- c k4. c q4. c q'4.i C k't)) iw " 

+ jy 5Z((^ C H C q|. C q'ti C k't(^4) 

qq' 


In this paper, we adopt the decoupling procedure for 
Green’s functions including several operators (2^. I 25 I - I 27 I ] . 
For example, 


((<5+Ckt C q4. c q' t ; C k't)) - ( S '+ C qf C ‘l'tK( C kt; C lc't)) 

- ( 5 '+ C qt C kt)(( C q't; C k't))- ( 25 ) 

We need further the Green’s functions ((S^Ckj.; c^,^}} and 

((5_Ckf; 4, t )). We neglect the terms of the order of J in 
the equations of motion for these Green’s functions; this 
means that we use the following approximation: 


I/O t \\ 'Uxfo'X d* XVyhy .. t \\ 

((SzCk4.; c k' t ))iw n - ---7—— (W*Ckt;4't))"„, 

1 Vz^z ~r /Z 


(26) 




((•S-Ckj,; c k ,*))j 


(27) 


For F k k' = ((ck 4 ,i c k'|))ia;„, we neglect the terms of the 
order of J in a similar way to obtain 

r’’ (■ \ T i'Uyky) ^ /. \ /'oo\ 

Fkk'(* w ") — ” j ] j Gkk't'(* w n) T 0 (d). ( 28 ) 

IU) n + V z k z + H 


to k. Here we set 


m k 


= 3 Z]Xqt c kt5'-), 


Tiktr — 5 X C c|CT- C kcr)) 

q 


and define 

g£(*w„) = - 


iw n + /j 


(iw n + ^) 2 - (v 2 fc 2 + u 2 fc 2 + v 2 k 2 ) 
iu) n + n + av z k z 


° ka (Wn) (iw„ + m ) 2 - (« 2 fc 2 + ^ 2 + « 

We also define the following functions: 

F{iu n ) = ^E G k(* w ™)’ 

k 

FOn) = -^ E ( mk - G k(®n), 


(30) 

(31) 

, (32) 

• (33) 

(34) 

(35) 

(36) 


G(iu n ) = — 55 ( n kt + «kt - 1) Gl(iu n ). (37) 


Then we obtain 
J 


E r kk'T = J 


r (iio n )G^t^.(iuj n ) 


2 1 + JG(ico n ) + (^) r (iu n )F(iu n ) 


(38) 


Now we can obtain a closed solution for a set of Green’s 
functions. The Green’s function Gkk'tX^n) reads 

-<o (■ , \ J rio / ■ ' J ^ 


Gkk / '|'(^ri) — f ^kk , Gk^(^CU n ) 2y Gk|(^) 2 ^(^ n ) 

x Gk^Xcu n )- 


1 + JG(iijj n ) + (ij) r(jw n )F'(iw ra ) 


(39) 

The Green’s function Gkk't(* w n) is similarly obtained as 


We set (S z ) = 0. From eas. El) and (El) . Tkk't is written 
as 


Gkk't(^n) — ^kk'GkX^n) 


2^G^(iWn)fr(iw n ) 


r kk'f(^n) 


iu n + n 


(iu n + n) 2 - (v 2 x k 2 x + vlkl + v 2 z k 2 ) 


m k ) r> »T G k't(® W ") 


- m k 


27V k ' f 
2 TV 






- (n k t + nkj. - l) ^ E r pk't(* w «) 


(29) 


x Gk/i(*w n )-j-. 

1 + JG(iuj n ) + (-j) r(zw n )F , (iu; ra ) 

(40) 


IV. KONDO TEMPERATURE 

) 

From the Green’s function Gkk"|-(*w„), the Kondo tem¬ 
perature 7k is determined from a zero of the denomina¬ 
tor in this formula. We perform the analytic continuation 
iui n —» w and consider 


where • • • indicates the terms, being proportional to k z . 
which give small contributions to G k k't because they 
would vanish when we take the summation with respect 


1 + JG(w) = 0 (41) 

in the limit ui —> 0 by neglecting higher-order term being 
proportional to (J/2) 2 . We neglect the term of the order 
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of J in Tikcr = (cocker)- The equation in ea. fllT) is written 
as 


1 = 


!>£<“> 

k 


tanh 


( e k ~ M 
V2fcBTk 


+tanh 


-fk ~ M 

2fc B X K 
(42) 


where £k = y u 2 fc 2 + u 2 fc 2 + v 2 fcf. Let us adopt for sim¬ 
plicity that v x = v y = v z = v and the dispersion is 
£k = ±vJ fc 2 + ky + fc 2 — p,. We can consider, in gen¬ 
eral, the d-dimensional case; by setting v z = 0 for two 
dimensions and v y = v z = 0 for one dimension. Then 
the equation for Xk is 


1 = I J-i£ f dkk d ~ x f---- 

8 (2n) d J 0 \uj + p — vk 

1 

K \vk — u — i7r(2n + 1 )Xk 

1 

u/c + fj, + iir(2n + 1 )Xk -I ’ 


w + /x + vfc 


(43) 


where fid is the solid angle in d-dimensional space, 
namely, the area of the (d — l)-sphere S'* -1 . D is an 
cutoff and we set fc B = 1 for simplicity. When d is odd, 
the integrand is an even function of k. In this case the 
integral with respect to k is reduced to an integral in 
the complex plane. For d = 3, the equation in eq. IS is 
written as 


1 = 


32tt 2 


J 


[E; 

n> 0 


1 ( p + i-K(2n + l)Tk \' 


x 27 tiTk 


1 


1 


-£ 

n< 0 

x 27 tzT k 


lo — in(2n + 1 )Tk to + 2p + iir(2n + l)Tk 
1 f /i + iir(2n + 1 )Xk n 2 
v 


1 


w — *7r(2n + l)Tk ox + 2/x + iir(2n + 1)Tk 

(44) 


where the summation has the upper limit no = 
D/(2ttTk)- By using the formula of the digamma func¬ 
tion, 


n o 

£: 

n=0 


= ^ o+a; + ^0 -V , o+ :r K ( 45 ) 


we obtain 


1 = 


32tt 2 


1^1 




v 1 2 2 tHT k 


+ n 0 - V’ ( o - 


27tzTk 


,1 cj + 2/i \ /1 cj + 2/i\ 

_ ^ 2 + 27riT K + n °) + ^ ( 2 + 27riT K ) 


, . 1 w 
+ ^ 2 + 27riT K 

1 2 27txTk 




n 0 


.1 cj + 2/i 

[2 ~ 2niT K 


(46) 


We assume that the cutoff £> is much larger than the 
temperature: D ^$> T. We employ the asymptotic form 
of ip(z) ~ log(z) for large 2 Then the Xk is the solution 
of the equation 


1 = jj;Pd(p)\ J \ 


AT(f) +log 


D 


^W+Ap? 


where we defined t = T/\p\ and 

We introduced the density of states pd as 
(27T) a \ V J V 


, (47) 


(48) 


(49) 


Here, fid is the solid angle in d-dimensional space, that 
is, the area of the {d — l)-sphere S d ~ l . As t approaches 
0, K(t) behaves as K(t) ~ log(l/ 7 rt). For large t, K(t) ~ 
7C(3)/(t r 2 f 2 ) where £(3) is the Riemann zeta function at 
argument 3. The equation eq. dTTl) always has a solution 
for |J| > 0. When pd|T| is small, Pd|T| <C 1, we have 
a solution in the logarithmic region of K(t). Then we 
obtain the Kondo temperature, 


& B 7k 


2e 7 X> 1 



exp 


8 ^ 
p^m)' 


(50) 


where 7 is Euler’s constant and fc B is included in the 
j formula of Xk- The result shows that the Kondo effect 
indeed occurs in a Dirac system. 

When D |/z|, Xk is 


fc B X K = 


4e 7 |/x| 
7 r 


exp 


-j_). 

pd(n)\j\J 


(51) 


In the limit |/i| —>■ 0, the equation 1 + JG(0) = 0 has 
no solution. This indicates that when \p\ is small, Xk is 
reduced and vanishes for \p\ —> 0. Hence, when the Fermi 
surface is point like, the Kondo effect never appears; this 
is because the scattering from the localized spin becomes 
weak for the point-like Fermi surface. This is consistent 












































5 


with calculation obtained by using the Abrikosov-fermion 
mean field theory for a topological insulator j2a_ and that 
by the functional-integral saddle-point theory [29(. 

We also examine the case where p&\J\ is large being 
of order 1 although our method is likely, however, not 
reliable in this region. In this region, we have a solution 
for K(t ) ~ 7£(3)/(7r 2 t 2 )-|-0(l/t 4 ) and Tk is an algebraic 
function p such as p a with a constant a. If we use K(t ) ~ 
7C(3)/(t r 2 t 2 ), we obtain 




i.rmAja 

7 T 2 V 16 v\v 


(52) 


Tk is proportional to \p\ times an algebraic function of 

Pd\J\- 

In one dimension (d = 1), Tk is obtained in a similar 
way. For small p d |J|, we have 


&bTk 


2 e<D 


1 + 


D 2 


: exp 



For small p, Tk is 


&bTk = 


4e 7 |/i| 

7T 


exp 



(53) 


(54) 


Let us turn to the two-dimensional case. The integrand 
is an odd function of k for d = 2, and thus the integral 
is not straightforwardly reduced to a complex integral. 
The results for d = 3 and 1, however, show that fc d_1 
in the integrand is approximately replaced by {\p\/v) d ~ l 
because the zero of denominators in the integrand gives 
important contributions. This results in the formula of 
T k for d = 2. 

As a result, the formula of Tk in d dimensions reads as 
in ea. rtbUll for d= 1, 2, • • ■. The solution for Pd\J\ being 
of order 1 in the region K(t) ~ 7£(3)/(7r 2 t 2 ) is given by 
ea. (l52l) . 

We derived Tk by using the Green’s function theory. 
Tk is proportional to \p\ with the exponential factor for 
small pd\J\ and \p\ <C D. The formula of Tk shows that 
Tk decreases as the dimension d is increased. We sum¬ 
marize the results for Tk in Table 1 (the first column). 

There has been a calculation based on the numerical 
renormalization group method [30j for a pseudogap U = 
oo Anderson model with the density of states p(w) oc 
|oj + p\ r for r = l[Hl. Their results indicate a particle- 
hole asymmetry showing Tk oc p x with x os 2.6 for p > 0 
and Tk oc \p\ for p < 0. This kind of asymmetry cannot 
be visible in our method. 


V. ELECTRICAL RESISTIVITY 


We consider the conductivity given by the formula: 
2 e 2 


cr = — - 


' 2 d f 

TkV k -Q — p d d£k- 


(55) 


TABLE I: Characteristic behaviors of Kondo temperature 
Tk, the resistivity R and the specific heat AC where D = 
T/V1 + D 2 /(2p) 2 . We set K(t) = Re^(l/2 + ip/(nT)) - 
ip) 1/2) and I(t) = lrmp'(l/2 + ip/(irT)), where ip(z) is the 
digamma function, a is a constant. We assume that T<D. 


Quantities 

Conditions 



pd\J\ < 1 

Tk — D exp ( 

T k 

(M « D) 

T K ~Hexp( m( ® )|i7| ) 


Pd|J|~0(l) 

Tk — \p\(y/ pd\J\ + • • •) 


Tk < T <§; |p| 

(loe(i)) 

R 


— (l + |pa|d| log y ^ 


T ~ \p\ 

R^^\J\ 2 {l + ±Pd\J\K(t)) 


\p\ « T 



Tk < T <§; |p| 

AC<x(ip d |J|) 3/2 r(log^)' 3/2 

AC 


«(|pd|J|) 3 T(l + Ap d |J|logf) 


T~H 

AC oc T-§^ (CiTA'(t) + C 2 T/(t)) 


H « T 

AC ~ Q(p/T) 


The life time r k is given as 

— = cA^ImGkkt(w + id) -1 , (56) 

Tk 

where c is the concentration of magnetic impurities. We 
adopt that T > Tk so that mk = 0. Then we have 
r(w) = —(3/4 )F(u>) and the Green function is 


GkktM = G° t H 




(57) 


1 + JG(u)-(j) |F( W ) 2 
This is approximately given as 

(58) 


ImF(w + id) is the density of states: 

ImF(w + id) = — -7rpd(w + p), 


(59) 






























6 


where we assume that p> 0. Then the life time Tk reads 

32 


Tfc(w) ~ 


3c7T J 2 p d (oj + fJ.) 
This results in the conductivity: 


(1 + JG(u)). (60) 


2e 2 2 32 / 1 ... _. 

-v - -—I l + —pd(^)\J\ 


3cnJ 2 


-2i/j 


D 

2nT 


+ 2ip 


'1 p D 

ib - + -1- 

v 1 2 iff 2 t tT 


-<H 2 


— ) 

7 xiT J 


■iP 


7 xiT 27 tT 




niT J 


(61) 


We assume T <C D. At low temperatures, when T <C 
11 , we obtain 


This has no log-T dependence. As a result log-T appears 
in the region Tk < T -C \p\. We summarize the results 
for R in the second column of Table 1. The log-T depen¬ 
dence is reduced as the chemical potential |/z| approaches 
the Dirac point, that is, log-T will disappear as p —» 0. 


VI. SPECIFIC HEAT 

Let us examine the specific heat in this section. The 
specific heat also shows a singularity as in other physical 
quantities. We calculate the additional entropy coming 
from the s-d interaction with magnetic impurities. The 
expectation value of the interaction part H s d is given by 


V = (H sd ) 

= 2 7V/3 ^ rkkV (* Wn ) 

kk'ncr 

J _ o J 

= — 5 GlAiUnWun) = Gi(iu n )t(iu) n ), 

k ncr kn 

( 66 ) 


1 


lo g ( jr) ■ ( 62 ) 


9c7T |J| 
The resistivity R = 1/tr is 


T \ 


R = Tl ci: ~F^\ J \ 2 
64 e z v z 




64 e 2 v 2 


1 + «(|H log (isT) 


■ (63) 


with D = D/yJ 1 + D 2 /(2p) 2 . There is a logarithmic 
T-dependence in the resistivity, which characterizes the 
Kondo effect. This is consistent with the result in Ref. [32] 
up to log-T term. We must note that the coefficient of 
the logarithmic term depends on the chemical potential 
M- 

When T ~ 0(\p\), R reads 


i? -^ C7r ^ |J|2 ( 1 + ^ d(M )l J l 


K(t) 


T log 


D 


y/D 2 + 4 p 2 


(64) 


The log-T dependence is reduced in this region. In the 
high-temperature region where |/x| <C T holds, the log(T) 
term will be absent. In the high-temperature region for 
\p\ <T, the resistivity is 


1 + ywMi-'i 


7 <(3)^(f) 2 -3i«5)T( 


T 


I- (65) 


where t(z) is defined as 


J 


T( 2 ) 


'‘ : ' ~ 4 1 + JG{z) + {J/2) 2 T{z)F{z) ‘ 

The a) n -summation is performed to give 


(67) 


F = 2jly4- 

v 27T* 


Re 


dujf(w) 


— id) 


uj + p + vk u j + /j, — vk 

+ i7T^S(uj + i-i + vk) + d(u> + p — vk)^Ret(ui — id) 


( 68 ) 


Because the w-dependence in t(w) is important, we ne¬ 
glect the w-dependence in l/(w + p ± vk) and d(ui + 
p ± vk). We adopt to*, = 0 (T > Tk), so that 
r(w — id) = —(3/4 )F(oj — iS). F(co — id) is approxi¬ 
mated as F{ui — id) ~ F{—id) ~ pd{p)F d + inpd(p)/2, 
where F d = — sign(/z)i>A/|/i| for d = 3 and F d = 
— sign(/r)(l/2) In |(v 2 A 2 — p 2 )/p 2 \ for d — 2 with a cutoff 
A in fc-integration. The real part of F(u> — id) is propor¬ 
tional to p/\p\. Then the interaction energy is 

v = + < 69 > 


Using the relation between the free energy and V given 
as 


V = J 


dF 

dJ’ 


(70) 
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the additional entropy is obtained as 


d 


J dJ' d 


AS = -—(F-F 0 ) = - I ——V(J',T). (71) 


J' dT 


We employ the expansion formula for the Fermi distri¬ 
bution function /(w): 

f D n 2 

/ dujf(oj)h(uj)= / dujh(u}) +—(kBT) 2 h'(0), (72) 
Jd Jd 6 

for a differentiable function h(w). Then we have 


V ^-F^p^j) 2 J 


du- 


1 


6 (fcsT)2 &cl + JG(w) 


—£) 1 + J G(cj) 

1 


cj—0 


(73) 


We first examine the low-temperature region character¬ 
ized by T -C \g\. When \g\/T is large, the second term 
of V is written as 


V 2 = - Y F dPd (f,)\J\(k B T) 2 


1 


(ln(T K /T)- 5 (o;)) 




where g(ui) is defined as 
sM = In (+ 


8 


TJ P d(ji)\J\ 

We have g( 0) = 0 and g'(0) is evaluated as 


cj—0 

(74) 


(l + JG(w)) ■ (75) 


ff'( 0) ^ - 


1 


2niT 


u + 




i 


2niT 
1 

2iriT 

1 


1 

2 niT J 2iviT 

D 




1 _ k_\ 

2 niT J 


$ ( - - — , 

1 2 niT 2nT 

I 1 M -D 

1 2 7 t?T 2nT 


D 2 


2p D 2 + Ap ?' 


(76) 


This is in contrast to the case of conventional Kondo 
effect where </(0) is proportional to the inverse of tem¬ 
perature 1/T. Then we obtain 


V 2 = ^A D p d (g)\J\^(k B T) 2 
4 H 


(HTk/T)) 


2 ’ 


(77) 


with Ad is a constant given by D 2 /(D 2 +4 p 2 ). The first 
term of V has a similar singularity[22|. We expand this 
term in terms of the inverse of ln(Tk/T) to obtain 


Vi = -3F dPd (n)\J\ 


1 


(ln(T K /T)r J-d 


du}g(oj) 


( 78 ) 


g(oj) is expanded by means of ui/T by restrict¬ 
ing the integral region to (— k B T, 0) to give a term 
(k B T) 2 / (ln(Tk/T)) 2 . This contribution is included in the 
coefficient Ad such as Ad = (1 — 3/n 2 )D 2 /(D 2 + Ag 2 ) 
and we have 

U = ^A DPd (g)\J\^(k B T) 2 - -1--j. (79) 

M (ln(T K /T)j 


Because of the relation lnTk = ln(2e 7 -D/7r) + 8/(p d J), 
the following holds: din 7 k = —8d(pdJ)/(pdJ) 2 - Using 
this, the contribution to the free energy is evaluated as 


A F = 


rPdJ 


d(pdJ) 

PdJ 


V(pdJ) 


2 tt 2 A d ^ (k B T) 2 ^ 5 / T ))2 


PdJ 


1 


ln(T K /T) ln(D/T) 


In 


Pd J ln (TK 
T 


, (80) 


where we neglected the factor 2e 7 /7r for simplicity. The 
entropy AS = —<9A F/dT is found to be 


AS = An 2 A D k B F d — 
A 4 


1 

3 



1 

2 



3 fp d J\ 

5 l 8 J 





(81) 


The entropy has a singularity such as T In T in the fourth 
order of pdJ , which is rather weak singularity compared 
to the conventional Kondo effect with InT term. This 
results in the specific heat AC = TdAS/dT: 


AC „k B Tfp d J\ 

- — AoFi — (t) 




2 V 8 


T 


3/2 


(ln(T/T K )) 


3/2 ’ 


(82) 


Hence, the specific heat in a Dirac system has a singular¬ 
ity as a function of the temperature near Tk- This should 
be compared with the specific heat anomaly in the orig¬ 
inal Kondo effect given as AC/kn ~ l/^ln(Tk/T)^ , 
which is shown by following the above method for the 
original s-d model. The singularity in a Dirac system is 
much weaker than that of the original Kondo effect. The 
log-T appears in the range Tk < T <C \g\/k B . 

Let us then examine the specific heat in the intermedi¬ 
ate region T ~ \g\. The expectation value of the interac¬ 
tion term V in ea. (173l) is evaluated by expanding JG(lo) 







































iii w and the integral is restricted to the interval (—T, 0). 
JG(lo) is written as 

JG(w) = JG( 0) + ^Pd|Jl^Inw/'' Q + i^j , (83) 

and the main contributions to V are 


V~F d ( Pd \J\) 3 4 5 6 7 


C\TK{t) + CzTImip' 


nT 


(84) 

where C\ and C 2 are constants. This results in the spe¬ 
cific heat given as 


with the characteristic temperature scale 


l B T K = 5exp(--((W), (87 ’ 

for small p d (p)\J\ "C 1- This is the conventional form, 
as in the original Kondo problem, being proportional to 
exp(—const/p| J|) with the density of states p. When \p\ 
is small compared to the cutoff D , Tk is proportional to 





Pd(p)\J\J 


( 88 ) 


AG 1 B 2 

— ~ --F d (p d (p)\J\fk B T— (CiTK(t) + C 2 TI(t )), 

(85) 

where we defined 


In the range Tk < T -C \p\ /fcg, the log-T anomaly ap¬ 
pears in the physical quantities. For large p d \ J\ of order 
1, Tk will be given by an algebraic function Q{x ) such as 


I(t) =lmtf)' Q + i 1 ^fj ■ ( 86 ) 

Thus the log-T terms do not show up in the region T ~ 

H 

In the high-temperature region defined by |/x| -C T, the 
temperature dependence of V comes from the terms of 
order p/T. Thus the additional entropy is also of order 
p/T and is negligible in the region |/z| <CT. 

VII. SUMMARY 

We investigated the Kondo problem with the localized 
spin which couples to Dirac fermions. The Kondo tem¬ 
perature Tk was calculated from a singularity of Green’s 
functions. The logarithmic terms in the resistivity and 
the specific heat were derived in the low-temperature 
region, where The region fceT -C \p\ is called the 
low-temperature region. We considered two regions: 

(1) When /cbT -C \p\, the Kondo screening occurs 


fc B Tk — \p\Q(pd(p)\J\)- (89) 

(2) When fcsT is not much smaller than |/i|, namely, 
ksT ~ \p\ or fcsT > \p\, the Kondo effect is suppressed 
and the resistivity and specific heat do not exhibit a loga¬ 
rithmic (log-T) anomaly. Tk vanishes when the chemical 
potential is just at the Dirac point, leading to the absence 
of Kondo screening. 

The vanishing of Tk when p is at the Dirac point 
is consistent with the results for the pseudogap Kondo 
problem [Tsj, and suggests the existence of a phase tran¬ 
sition, as |J| is increased, from non-screening phase to 
screening phase. The term being proportional to a z like 
the magnetic field would also be important in the Kondo 
effect, which we did not consider in this paper. It is 
also interesting to study the nature of interaction be¬ 
tween two magnetic impurities in Dirac metals. This 
issue was investigated intensively in the original Kondo 
problem [3^-0]. These are future subjects. 

The author expresses his sincere thanks S.A. Jafari, I. 
Hase and K. Yamaji for discussions. He also thanks M. 
Vojta, E. Orignac and A. Principi for useful comments. 
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